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Abstract 

The hyperfine structure of hsnd 1 £>2 states in 87 Sr was 
studied for principal quantum numbers 11 < n < 25 
by Doppler-free two-photon excitation combined with 
thermionic detection. Based on the more accurate 
data, compared to previous results, the hyperfine 
structure was analysed by second order perturbation 
theory using MQDT wavefunctions. As a result we 
find that the surprisingly large B-factors obtained by 
fitting the data with the Casimir formula do not re- 
flect a quadrupolar interaction but rather consist pre- 
dominantly of second order magnetic dipole contri- 
butions. Further, inconsistencies in the King plot of 
odd-even isotope shifts were removed by this type of 
analysis. The result demonstrates the importance of 
higher order effects in hyperfine splittings of excited 
states in atoms with two valence electrons, even when 
the fine structure splitting is large compared to the 
hyperfine splitting. 

PACS numbers: 35.10.Fk, 32.80.Rm, 31.30.Gs 



1 Introduction 

Many problems in atomic physics can be treated 
using perturbation theory. Normally for low lying 
states with large energy distances to neighbouring 
levels the application of perturbation theory in first 
order gives satisfying results. However, the level spac- 
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ing of high lying states (Rydberg states) is consider- 
ably smaller and higher order perturbation theory is 
required for correct treatment. 

The electron configuration of excited states of the 
alkaline earth elements is determined by the two elec- 
trons outside a core of closed electron shells. The 
existence of low lying bound doubly excited states 
causes configuration mixing and influences the cou- 
pling of the two valence electrons (singlet-triplet 
mixing). The multichannel quantum defect theory 
(MQDT) has proved to be an efficient tool jlj for the 
description of such perturbed Rydberg states with 
only a small set of parameters. 

Particularly the atomic hyperfine structure (hfs) is 
sensitive to admixtures in the wave function of Ryd- 
berg states Q and is therefore a most sensitive quan- 
tity for investigating such phenomena. 

For the evaluation of the hfs it is in most cases suf- 
ficient to consider only the magnetic dipole and the 
electric quadrupole interaction. Then it is possible to 
parameterize the hfs - splitting in first order pertur- 
bation theory by the hyperfine constants A and B, us- 
ing the Casimir formula. The hyperfine constants A 
and B measure the magnitude of the magnetic dipole 
and the electric quadrupole interaction, respectively. 
This kind of evaluation and separation of the interac- 
tions assumes that contributions of higher order per- 
turbation theory are negligible. If that is not the case 
the hyperfine constants lose their physical meaning. 

Previously, the effects of higher order contributions 
to the hfs in Rydberg states due to neighbouring 
fine structure levels have been observed in Hg 
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3 He jpl, Ca §, SrJ,0-0, Ba Sm @, 

Yb [18|, and Pb |l7|. For very high Rydberg states 
the energy spacing of states with different principal 
quantum number n is small enough to cause hyper- 
fine induced mixing [[l9|-pl"|. 

In the present work we report on high resolu- 
tion laser spectroscopic measurements of the hfs of 
5snd 1 £>2 states in 87 Sr and of transition isotope 
shifts of the stable isotopes 84 > 86 ' 88 Sr for principal 
quantum numbers 11 < n < 25. These measure- 
ments were done with enhanced accuracy compared 
to previous results 



12 , 23 , in order to obtain suffi- 
ciently precise data for accurate comparison of the 
hfs with the results derived from second order per- 
turbation theory. The hfs was calculated by means of 
MQDT-wavefunctions in first and second order per- 
turbation theory. We found that second order contri- 
butions need to be considered in the analysis of the 
hfs in order to explain the unexpected large hyperfine 
constant B and the shift of the hyperfine multiplet, 
which is responsible for large deviations from the lin- 
ear dependence in a King Plot analysis. 
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2 Experiment 

The experimental set-up used in the present work 
has been described in detail in an earlier paper pjj ]. 
Briefly, it consisted of a cw narrow bandwidth ring 
dye laser (Spectra Physics Model 380D) pumped by 
an argon ion laser (Spectra Physics Model 2045). 
Typical operating conditions were uv pump powers 
of 2.5-3.0 W, resulting with a stilben 3 dye in 200- 
400 mW single-mode power. The dye laser is tunable 
between 410 and 470 nm. With these wavelengths 
Sr I Rydberg levels above n — 8 can be excited from 
the 5s 2 ground state by Doppler-free two-photon 
excitation. A Michelson-type wave meter (Burleigh 
WA 20) with 300 MHz accuracy was used for iden- 
tifying Rydberg levels up to n = 200. For measur- 
ing the hfs - splittings the laser frequency was sta- 
bilized to an external confocal cavity. This, in turn, 
was controlled by an active molecular frequency stan- 
dard with an accuracy of better than 100 kHz, using 
several precisely known 1% levels combined with the 
frequcncy-offsct-locking technique |25|. A computer 



Figure 1: Two-photon excitation spectrum of the 
5s 2 — > 5sl2d 1 L>2 transition as a function of the 
energy detuning relative to the signal of the isotope 
88 Sr. 

(HP 21 MX-E) was employed for real-time data ac- 
quisition and for scanning and controlling the laser 
frequency. 

The Sr-atoms were excited and, following colli- 
sional ionization, detected with high efficiency in a 
thermionic ring diode (26). In order to improve the 
signal to noise ratio, the retroreflected laser beam was 
chopped at a frequency of 80 Hz for lock-in detection 
(EG&G PAR Model 5210). 

A typical spectrum of the 5sl2d 1 £>2 state is shown 
in Fig. 0. To determine the line positions, Lorentzian 
line profiles were fitted to the experimental data. 
The accuracy of the line distances was enhanced by 
recording and evaluating several spectra (up to 40) 
for each state. The mean values of the transition fre- 
quency differences with regard to 88 Sr are listed in 
Table |. 
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3 Theory of hyperfine struc- 
ture 

The hyperfine interaction describes the interaction 
between the electrons and the magnetic and electric 
multipole moments of the nucleus, where the even 
magnetic and the odd electric multipole moments 
vanish. In most cases it is sufficient to take into ac- 
count only the interaction with the nuclear magnetic 
dipole moment fij — <?/I and the electric quadrupole 
moment Q. According to Sobel'man j27|, the Hamil- 
tonian can than be expressed as 



H his = H » + Hc - 



(1) 



where 

^ = E^( £(1) -^o[ s(1)xC(2) )] (1) ) 

i 

^Q--e 2 E^" 3c(2) 'Q (2) ' ( 3 ) 



(2) 



and v- 1 ' and are the single electron angular mo- 
mentum and spin operators, respectively. The tensor 
C^ 2 ) is connected to the spherical harmonics Yi m by 
Ci = [M(2£+l)] 1 /^ m . 

In first order perturbation theory the hyperfine in- 
teraction yields an energy splitting 



A4 1 ' = ( 7 JIF\H Ms \jJIF) 



(4) 



of a state with total angular momentum F = I + J. 
Here I and J are the angular momenta of the nu- 
cleus and the electrons, respectively. Evaluating this 
matrix element, one obtains the well known Casimir 
formula $M 



(i)_l A1 s l D zK(K + l)-J(J+ 1)7(7+1) 



AE))' = -AK + B 
F 2 



27(27- 1)J(2J- 1) 



where 



K = F(F + 1) - J( J + 1) - 7(7 + 1) 



(5) 



(6) 



The hfs constants A and B describe the magnitudes 
of the magnetic dipole and electric quadrupole inter- 
actions, respectively. Due to the different dependence 
on F of the two terms in Eq. (||), it is possible to fit 
the Casimir formula to experimentally determined hfs 
- splittings, and to separate the contributions of the 
magnetic dipole and the electric quadrupole interac- 
tion. 

The measurement of the isotope shift (IS) is also 
possible since for first order perturbation theory the 
center of gravity of a hyperfine multiplct remains un- 
shifted E P AE^(2F + 1) = 0). 

The situation changes if one considers contribu- 
tions of second order perturbation theory. The sec- 
ond order energy splitting of a state \jJIF) is 



a4 2) = 



E 



|(7J7F|i7 hfs | 7 'J'7F)| 5 



{Ejj — E. 



i'J') 



(7) 



Hence when second order effects are not negligible, 
i.e. when neighbouring fine structure levels with equal 
quantum number F are lying close together, the en- 
ergy splitting can no longer be approximated by the 
Casimir formula. Then it is not possible to sepa- 
rate magnetic dipole and electric quadrupole contri- 
butions from the experimental hfs - splittings, be- 
cause the second order contribution of the magnetic 
dipole interaction has a similar dependence on F as 
the first order contribution of the electric quadrupole 
interaction. This is especially important when the 
hfs is dominated by the magnetic dipole interaction. 

The second order energy splitting goes together 
with a shift of the center of gravity, since the con- 
tribution of each fine structure level has a different F 
dependence (only states with equal F are connected) . 
This has the consequence that the IS of a perturbed 
hyperfine multiplet can only be extracted using sec- 
ond order perturbation theory. 

To study the influence of the second order contribu- 
tions, one can calculate the hfs in first and second or- 
der perturbation theory if the fine structure splitting 
is known. For a two electron atom like Strontium the 
magnetic dipole and electric quadrupole interaction 
can be written in terms of the one electron matrix 
elements {ap) and {be). This has been worked out 
by Lurio et al. |29] for a s^-configuration. Explicit 
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expressions for a more general ^i£ 2 -configuration are 
given in the Appendix of the present work. The one 
electron matrix elements depend on the expectation 
value of r -3 for i ^ and on the probability density 
at the nucleus, |-0(O)| 2 , for s-electrons. According to 
Sobel'man j27| (ai) and (bg) can be written in units 
of 1 Rydberg (Roo) as : 



(ae=o) 
(bt=o) 



= gia 2 — {R n ^\R nl )R 
m p r 6 



2— o (RrU I -| \RrU)Roo 

^a 2 9l —al\m\ 2 Ro 







(8) 

(9) 
(10) 
(11) 



Here gi is the nuclear g-factor, a the fine struc- 
ture constant, m/m p the ratio of electron and proton 
mass, e the elementary charge, R n £ the radial part 
of the one-electron wavefunction and do the Bohr 
radius. To a first approximation, (R n e \o-Q/r 3 \R n i) 
and (0)| 2 can be calculated using the formula from 
Ref. Ml : 



(R n e\^\Rne) 



1 



£{£■ 



1 



i z 



Roo 



na'A rr 



3/2 



3/2 



(12) 
(13) 



where ?i* 3 is the effective principal quantum number, 
e n i the binding energy (in Rydbergs) of the nl orbital, 
and Z the effective nuclear charge. 

4 Calculation of the hyperfine 
splitting with MQDT wave- 
functions 

As shown in reference the MQDT is a powerful 
tool to describe perturbed Rydberg series of alkaline- 
earth atoms. 



The empirical MQDT study by Esherick [|o| 
showed that the J = 2 bound-state spectrum of 
Strontium can be described by a five channel MQDT 
model which includes the recoupling of 1 D 2 and 3 D 2 
channels and the perturbation of the 5snd 1 D 2 and 
3 D 2 series by the Adn's l D 2 , 4dn's 3 D 2 and 5pn"p 1 Z? 2 
channels. Here, n' and n" are the principal quantum 
numbers of the Rydberg series converging to the Ad 
and bp ionization limits. Within the accuracy of that 
analysis it was found that it is not necessary to in- 
clude the Adn"'d l D 2 channel. 

Since there are only a few doubly-excited states 
below the first ionization limit, the empirical MQDT 
study cannot give the exact admixture of each dou- 
bly excited channel to the levels although the total 
amount of the doubly-excited channels can be well 
distinguished from the singly-excited ones. In partic- 
ular Esherick's analysis can not be used to distinguish 
between the 5pn"p X D 2 and 4dn'"d 1 D 2 channels. On 
the other hand, MCHF calculations performed by As- 
pect et al. pl[ a nd ab initio MQDT calculations of 
Aymar et al. |32|] indicate that the Mn'"d component 
is twice as large as the 5pn"p one. Thus the wave- 
function of the even J = 2 bound states of Sr can be 
written as an expansion of pure LS states with the 
following terms 

| 7 J = 2) = Zl \5snd 1 L> 2 )ls + z 2 \5snd 3 D 2 ) LS 
+ z z \Adn's 1 L> 2 )ls + z 4 \4dn's 3 D 2 ) LS 

l ^\5pn"p 1 D 2 ) LS 

2 ^dn'"d l D 2 ) LS 

where the Zi are the mixing coefficients obtained by 
using Esherick's (3^] MQDT parameter. The mixing 
coefficients for the 5snd 1 D 2 and 5snd 3 D 2 states are 
drawn in Figs. ^| and ^ for principal quantum numbers 
11 < n < 25. The data show a strong singlet-triplet 
mixing in the 5snd configuration around n — 16. The 
main contribution of the doubly excited channels in 
this region of principal quantum numbers comes from 
the Adn' s configuration, where also a singlet-triplet 
mixing occurs. 



Z5 I 



As studied earlier 22 1 , the hfs of the 5snd 1 D 2 and 
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3 D 2 states is dominated by the Fermi-contact term 
of the magnetic dipole interaction of the 5s electron 
and depends strongly on the mixing coefficients z\ 
and 22- The n£ Rydberg electron has only a weak 
interaction with the nucleus. For pure hsnd Y D 2 Ry- 
dberg states the hfs - splitting should almost vanish 
and for pure hsnd 3 D 2 the dipole constant A should 
be nearly a^ s /\2. Only a weak quadrupole interac- 
tion due to the small admixtures of doubly excited 
channels is expected. 

In the calculation of the hyperfine splitting of the 
hsnd X D 2 states in second order pertubation theory 
states of different J and equal F enter the matrix 
elements (cf. Eq. [jj). These are the fine structure 
levels hsnd 3 Di, 3 D 2 , and 3 D% whereas the hsnd X D 2 
and 3 D2 states are described by Eq. |l4|, the expansion 
of the hsnd 3 Di and hsnd 3 D^ states are taken from 



the two channel MQDT analysis given in Ref. |34| in 
which the mixing of the hsnd and Adn' s configuration 
have been considered. Thus the wavefunctions of the 
J = 1 and J = 3 states can be written as: 

| 7 J= l) = a^snd^Ls + a 2 \4dn's 3 D 1 ) LS (14) 
\ 7 J = 3) = ^\hsnd 3 D 3 ) LS + p 2 \4dn's 3 D 3 ) LS (15) 

The mixing coefficients on and /3, can be calculated 
with the MQDT parameters given in Ref. |p3 . 

The evaluation of the matrix elements is done here 
without considering matrix elements between differ- 
ent configurations l\l 2 . This is a reliable procedure 
because these matrix elements vanish in most cases, 
especially for the Fermi-contact term, or are very 
small due to the different effective quantum numbers. 
However, in special situations Q the off-diagonal 
configuration matrix elements can be quit large. 

Further the energy differences between the 
hsnd X D 2 states and the other fine structure terms are 
needed for the calculations. These are listed in Table 
[j] and illustrated in Fig. [|. It can be seen that the 
energy differences of the three fine structure terms 
change their sign at different principal quantum num- 
bers. As studied earlier |nj], the small difference be- 
tween the hsnd X D 2 and hsnd 3 D% series at n = 19 
results in a large shift of the respective hyperfine mul- 
tiplets. 

In the calculation of the one-electron matrix ele- 
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Figure 2: Admixtures of the pure LS'-states according 
to Eq. H to the 5snd l D 2 Rydberg serie of Strontium 
calculated with the MQDT-parameters from Ref. [ |30[ 
as a function of the principal quantum number n. 
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Figure 3: Admixtures of the pure LS-st&tes according 
to Eq. [lj to the 5snd 3 D 2 Rydberg serie of Strontium 
calculated with the MQDT-parameters from Ref. |}0| 
as a function of the principal quantum number n. 



Table 1: Level energies of bound 5snd l D 2 , 3 D t , 3 D 2 
and 3 L> 3 states of Sr I (cf. Fig. 1 



E (cm" 1 ) 



E (cm" 1 ) 



E (cm" 1 ) E (cm" 1 ) 



11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 



44578 
44829 
45012 
45153 
45263 
45362 
45433 
45492 
45542 
45584 
45619 
45650 
45676 
45699 
45719 



6890 a 
6648 a 
0249 a 
2785 a 
6196 a 
1272 a 
2717 a 
6101 a 
2955 a 
1831 a 
7872 a 
2617 a 
5325 a 
3308 a 
2336 a 



44616.05 c 
44854.02 b 
45028.55 b 
45159.60 b 
45260.84 b 
45341. 36 b 
45414.43 b 
45475.35 b 
45527.10 b 
45570.97 b 
45608.37 b 
45640.43 b 
45668.11 b 
45692.08 b 
45712.94 b 



44620.08°" 
44860.06 b 
45036.95 b 
45171.49 b 
45276.65 b 
45350.35 b 
45420.84 b 
45479.88 b 
45530.18 b 
45573.28 b 
45610.07 b 
45641.68 b 
45669.06 b 
45692.90 b 
45713.50 b 



44625.1° 

44865.22 b 

45043.79 b 

45180.44 b 

45286.53 b 

45370.76 b 

45439.08 b 

45495.02 b 

45542.23 b 

45582.38 b 

45616.80 b 

45647.54 b 

45673. 10 b 

45695.94 b 

45715.80 b 



a Reference [ 
b Reference [ 
c Calculated 
d Reference [ 
e Reference [ 



value from reference 

53. 



merits, the experimental value of —1001(2) MHz from 
Ref. |39| was taken as O5 s -factor. The other one- 
electron matrix elements can be well approximated 
by introducing as s , and substituting Eqs. ( |l2|) and 
0) into Eqs. (|) to © : 



a 5« 



3/2 

a 5 



3/2 



£5s 



3 Q a 5 



1 



(16) 
(17) 

3/2 

(18) 



According to MQDT, the binding energies of the 
5s, 5p and Ad orbitals are assumed to be the same as 
those of the corresponding 5s, hp and 4d orbitals in 
Sr + . The additional binding energy of the nd, n's, 
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Principal Quantum Number n 

Figure 4: n dependence of the energy differences be- 
tween the 5snd 3 ^i,2,3 and the 5snd X D 2 states of 
Strontium. 



n"p or n"'d electron in a state can be approximately 
obtained by subtracting the energy of the inner orbit 
(eg. 5s) from the energy of the state \yJ). The ener- 
gies of the 5s, 5p and 4<i orbits in Sr + were taken from 
Rcf . |3^] , where the mean values of the fine structure 
components were used for the 5p and 4d orbits. For 
the quadrupole moment of 87 Sr, we used the value of 
0.335(20) • 10~ 28 m 2 from Ref. @, and for the dipole 
moment fi] — g]I = — 1. 089299 (l)/i at reported in 
Ref. The matrix elements in Eqs. (|l|) to @ 

were then introduced into the expressions listed in 
the Appendix to yield the energy splittings. 

The results of the calculations are shown in Fig. 
||. The second order contributions of the electric 
quadrupole interaction are less than 1 kHz and will 
be neglected here. It should be emphasized here that 
the first order electric quadrupole contributions are 
approximately one hundred times smaller than the 
second order magnetic dipole contributions. The sec- 
ond order contribution of the magnetic dipole inter- 
action in the 5sl9d 1 r>2 state is not shown here be- 
cause of the small energy separation to the 5sl9d 3 Da 
state (~ 2 GHz). In that particular case the contri- 
butions to the hfs are large in any higher order, which 
means that perturbation theory is not longer appli- 
cable. This is discussed in detail in Ref. ]To| . 

Due to the strong Fermi-contact interaction of the 
5s electron, the calculated total hfs - splittings for 
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Figure 5: Calculated energy-contributions to the hy- 
pcrfinc splitting of the 5snd 1 I?2 states in 87 Sr as a 
function of n. (a) : magnetic dipole interaction in 
first order, (b) : magnetic dipole interaction in sec- 
ond order (the very large contribution of the 5sl9d 
1 Z?2 is only indicated by the arrow and is discussed 
in detail in Ref. ]1C[]), and (c) : electric quadrupole 
interaction in first order perturbation theory. 
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these states depend almost entirely on the mixing 
coefficients Z\ and Z2 and on the experimental value 
of the o,5 S -factor. 



5 Comparison with the experi- 
mental data 

5.1 Hyperfine Splitting 

In Fig. H the calculated hfs - splittings are compared 
with the experimental data. The isotope shifts of the 
experimental data are taken into account by shifting 
the centers of gravity of the calculated hyperfine mul- 
tiplets to overlap with the experimental ones. Good 
agreement between experimental and theoretical data 
is achieved. The small deviations of the individual 
components can be explained by the fact, that the 
hfs is much more sensitive to the mixing parame- 
ters in the expansion of Eq. (ITj), especially to z\ 
and Z2, than to the energies of the states which are 
used in the empirical MQDT analysis (see Ref. p 7 ^]). 
To emphasize this point we have done a least square 
fit of the theoretical line positions for each state by 
varying the center of gravity of the multiplet and 
the mixing coefficients z\ and 22 under the condition 
z i + z 2 = constant. With only small changes of the 
mixing coefficients, which cause large effects on the 
hfs splitting, one can obtain a perfect agreement be- 
tween experimental and theoretical line positions, as 
is shown in Fig. The optimized mixing coefficients 
differ very little from the MQDT mixing coefficients 
except near n — 15, as can be seen in Fig. ^ and 
Table ^. Therefore, in the following discussion it is 
sufficient to use the MQDT mixing coefficients for the 
description of the main features of the hfs splitting. 

The need for considering second order effects for 
the hfs becomes evident, when one tries to separate 
the magnetic dipole and electric quadrupole contribu- 
tions by fitting the Casimir formula (which contains 
only first order effects) to the measured and calcu- 
lated hfs - splittings. This is shown in Fig. ||(a) and 
^(b) and is listed in Table |]. At the first glance an ex- 
cellent agreement between experiment and theory is 
obtained. However, the large constant B in Fig. ||(b) 
does not describe the electric quadrupole interaction, 
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Figure 6: Comparison between the measured and cal- 
culated (in second order) hyperfine splitting of the 
5snd 1 C>2 states in 87 Sr. The isotope shift is taken 
into account by shifting the centers of gravity of the 
calculated multiplets for agreement with the experi- 
mental ones. 
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and reflects mainly the second order contribution of 
the magnetic dipole interaction. The contribution of 
the quadrupole interaction is very small indeed, and 
is determined by the admixture of the Adn's config- 
uration. This is shown in Fig. |^(c), where the elec- 
tric quadrupole interaction is considered to first order 
only. 

5.2 Isotope Shift 

The influence of the second order effects on the IS can 
be seen in a King-plot analysis in Fig. [HJ The King- 
plot procedure (see Ref. f43|| ) can be used for testing 
the consistency of the IS measurement in two tran- 
sitions. It allows the separation of the specific mass 
and the field shift of a transition, provided the spe- 
cific mass and the field shift are known for a reference 
transition. In a King-plot, the modified IS £ A ' A of a 
transition i between two isotopes with mass numbers 
A and A' (where A' > A) is plotted against the modi- 
fied IS £ A ' A of the reference transition. The modified 
IS is defined as the difference between the experimen- 
tal IS Sv A ' A and the normal mass shift multiplied by 
the factor mArnA> / (rriA' — tua), where tua and ijia' 
are the masses of the isotopes : 



t A,A' r A. A' 



m A m A ' 
TUA' — m A 



m e Vi 



(19) 



Here m e is the electron mass and Vi the transition 
energy. In a King-plot the points for different isotope 
pairs should lie on a straight line. 

For the transitions 5s 2 ^So — > 5snd in Sr I 
this is true only for the even isotope pairs. As an 
example, Fig. |l^ shows a King-plot for the transi- 
tion Sr I 5s 2 l S -> 5sl7d 1 D 2 . The transition Sr II 
5s 2 S 1 /2 — > 5p 2 P\/2 from Ref. |39| was taken as ref- 
erence transition. Obviously the points involving the 
odd isotope 87 Sr (which are calculated using the cen- 
ter of gravity of the experimental hfs - splitting listed 
in Table P are far off the linear dependence of the 
even isotope pairs. This effect was also noticed in 
Ref. |1| for Ytterbium, where the deviation from a 
straight line was taken to determine the magnitude 
of the second order contributions. If one corrects the 
experimental IS of the 87 Sr for the second order con- 
tribution to the shift, represented by the center of 



8? Sr I 5s 2 1 S -> 5snd 1 D 2 

' ' : Experiment 

i 1 : Theory 











| F=5/2 7 / 2 


9/2 


11/2 


13/2 










F=13/2 


11/2 


9/2 


I 7/2 



i/ 87 -!/ 88 (MHz) 



Figure 7: Comparison between the measured and cal- 
culated hyperfme splitting of the 5snd x Di states in 
87 Sr. Compared to Fig. || the theoretical line posi- 
tions are optimized by a least square fit of the mix- 
ing coefficients Z\ , Z2 and the center of gravity of the 
multiplet. 
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Figure 8: Comparison between the mixing coefficients 
Z\ and Z2 of the 5snd D% states in Sr calculated with 
the MQDT parameters from Ref. |30| and optimized 
by least square fitting with the experimental hfs - 
splitting. The errorbars estimate the uncertainty in 
calculating the mixing coefficients from the MQDT 
parameters. 
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Figure 9: Hyperfine structure constants A and B for 
the 5snd 1 L>2 states in 87 Sr as a function of n (see 
Table ||). In (a) and (b): A and B are determined by 
fitting the Casimir formula to the measured and cal- 
culated line positions, (c): Theoretical prediction for 
B in first order. It should be noted that the constant 
B in (b) does not reflect the amount of the electric 
quadrupole interaction but rather the second order 
magnetic contribution. 
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gravity of the theoretical hfs - splitting listed in Ta- 
ble ^ (which is zero without second order effects) , the 
Table 2: Mixing coefficients zi and z 2 for the data fit a straight line with excellent agreement. 
5snd 1 £>2 Rydberg serie of Sr I (see Eq. Jl4|)) (a) 
calculated with the MQDT-parameters from Refer- 
ence p0[ and (b) optimized by least square fitting of 



the theoretical line positions (see also Fig. M. 



n 


{< 




(b) 


Z\ 




Z\ 




11 


0.971(3) 


-0.09(2) 


0.9732 


-0.0651 


12 


0.966(5) 


-0.14(3) 


0.9613 


-0.1658 


13 


0.956(2) 


-0.195(9) 


0.9540 


-0.2024 


14 


0.923(2) 


-0.309(6) 


0.9226 


-0.3111 


15 


0.86(3) 


-0.46(5) 


0.8252 


-0.5090 


1G 


0.81(4) 


0.54(7) 


0.7560 


0.6174 


17 


0.88(2) 


0.44(3) 


0.8771 


0.4400 


18 


0.943(6) 


0.29(2) 


0.9282 


0.3361 


20 


0.9667(9) 


0.220(4) 


0.9614 


0.2426 


25 


0.9847(2) 


0.152(1) 


0.9808 


0.1753 



? 0.3 

N 
X 

I— 

0.2 - 



• : uncorrected 
□ : corrected 
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(THz-u) 

Figure 10: King plot for the two-photon transition j: 
Sr I 5s 2 x 5 -> 5sl7d 1 D 2 . The values for the refer- 
ence transition i: Sr II 5s 2 S'i/2 — > 5p 2 Pi/2 are taken 
from Ref. |39|. The shifts for the isotope pairs in- 
volving 87 Sr are corrected by an amount given by the 
center of gravity of the theoretical hyperfine splitting 
listed in Table § 



6 Summary 



Higher order contributions to the hfs and to the IS 
have been studied for 87 Sr hsnd 1 Z?2 states in the 
range 11 < n < 25. An enhanced accuracy, com- 
pared to previous experiments, made it possible to 
identify second order magnetic dipole contributions 
to the hyperfine splitting and shift. 

In parallel the hfs was calculated in second order 
pertubation theory using MQDT wavefunctions. We 
demonstrate that the knowledgement of second or- 
der contributions is important for the determination 
of the quadrupole splitting and the IS, even when 
the fine structure separation is large compared to the 
hfs - splitting. Neglecting the second order contri- 
butions in the hfs analysis of the experimental data 
leads to unrealistically large quadrupole factors in the 
Casimir formula as well as to erroneous positions of 
the centers of gravity which do not give the correct 
IS. 

In the investigated region of n of the hsnd 1 D2 
states, it was shown that the second order mag- 
netic dipole contribution to the hfs - splitting is near 
n = 15 up to one hundred times larger than the first 
order electric quadrupole contribution which is deter- 
mined by small admixtures of doubly excited states. 
Because of the dominance of the second order mag- 
netic dipole contributions it was not possible to make 
any reliable predictions, based on experimental data, 
about the admixtures of doubly excited channels into 
the wavefunction of the hsnd x Di states that could 
significantly improve the results of a MQDT-analysis. 
Evaluating the IS from the 87 Sr spectrum, corrected 
for the second order contributions of the hfs to the 
shift, leads to a value consistent with a King-plot 
analysis. 
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8 Tables 



Table 3: Mean values of the frequency differences with regard to 88 Sr for the transitions 5s 2 ^So — > hsnd 1 £>2 
of the isotope pairs ( 87 Sr, 88 Sr), ( 86 Sr, 88 Sr) and ( 84 Sr, 88 Sr). 



n 






V s ' 1 - i> 88 (MHz) 






„m _ „ss (MHz) 


v si _ u ss (MHz) 


F = 5/2 


F = 7/2 


F = 9/2 


F = 11/2 


F = 13/2 


11 


185.83(45) 


89.81(25) 


-33.67(65) 


-179.95(28) 


-352.37(29) 


-244.99(21) 


-511.50(49) 


12 


634.105(95) 


395.214(82) 


88.69(12) 


-285.218(84) 


-724.297(87) 


-248.785(57) 


-521.06(13) 


13 


799.25(20) 


507.14(42) 


133.57(28) 


-320.94(25) 


-854.17(48) 


-240.73(11) 


-502.95(26) 


14 


1268.92(67) 


821.63(71) 


249.24(23) 


-439.12(24) 


-1240.48(72) 


-241.03(18) 


504.66(30) 


15 


2019.3(11) 


1372.1(14) 


523.41(80) 


-543.91(73) 




-246.54(83) 


-516.71(68) 


16 


-1844.34(65) 


-1299.56(56) 


-598.83(26) 


267.47(20) 


1295.54(66) 


-242.60(18) 


-508.06(30) 


17 


-1651.93(67) 


-1172.44(47) 


-550.56(21) 


217.64(27) 


1134.18(39) 


-234.36(14) 


-490.29(15) 


18 


-1404.38(42) 


-1019.48(50) 


-511.16(32) 


126.10(17) 


904.67(26) 


-230.05(15) 


-481.55(24) 


19 












-227.97(25) 


-477.05(42) 


20 


-1054.77(18) 


-712.84(23) 


-287.75(15) 


212.58(14) 


777.65(18) 


-226.23(14) 


-472.68(29) 


25 


-811.90(20) 


-566.69(23) 


-255.836(82) 


118.212(64) 


550.13(21) 


-222.179(49) 


-464.76(18) 
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Table 4: A and B constants and centers of gravity for the 5snd 1 i?2 states of 87 Sr determined by fitting the 
Casimir formula to the measured and calculated line positions (see Fig. ^(a) and ||(b). The theoretical line 
positions are calculated with the mixing coefficients obtained from the MQDT parameters given in Ref . [Q . 
It should be noted that the isotope shift is not taken into account for the theoretical values (Ai/|7). 



Experiment Theory 



n 


A (MHz) 


B (MHz) 


(i/g/ - O (MHz) 


A (MHz) 


B (MHz) 


(Ai/f/) (MHz) 


11 


-26.855 (20) 


4.49 (57) 


-111.72 (15) 


-36.06 


2.2 


-4.48 


12 


-67.895 (57) 


3.55 (15) 


-114.21 (10) 


-55.50 


3.2 


-6.11 


13 


-82.605 (26) 


6.59 (50) 


-112.43 (16) 


-79.39 


5.6 


-7.56 


14 


-125.217 (37) 


24.01 (79) 


-119.22 (24) 


-124.50 


23.1 


-13.96 


15 


-194.02 (18) 


-75.0 (32) 


-80.70 (99) 


-177.56 


-77.4 


22.74 


1G 


157.248 (34) 


14.85 (75) 


-121.80 (21) 


152.84 


15.2 


-16.36 


17 


139.542 (28) 


18.64 (61) 


-126.13 (17) 


139.57 


22.5 


-22.54 


18 


115.29 (13) 


45.96 (43) 


-150.29 (13) 


103.66 


48.2 


-48.68 


20 


91.067 (11) 


-50.62 (25) 


-29.50 (10) 


83.65 


-47.1 


71.22 


25 


67.954 (69) 


-17.12 (24) 


-56.78 (10) 


59.36 


-12.6 


47.38 
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9 Hyperfine interaction matrix elements 
9.1 Magnetic dipole interaction 



(e 1 £ 2 SLJIF\H^\£ 1 £ 2 S'L'J'IF) = (-1) J '+ /+F y/l(I + l)(27 +1) | ^ J' I } 

2 

x (^ 2 SI,J||^a^ (1) -\/To[ S « x C' 2 ')]' 1 ' +a Si ^ i ,o8 (1) ||^ 2 5'L'J') , 



i=l 



(^i^SL J| | ^ a u e {1) | li&S'L'J') 



{ _ 1)S +l'+j+ M ^ (2 J + 1)(2J' + 1)(2L + 1)(2L' + 1) <| ^ ^ j 



K)(-l) i V^i(^i + l)(2^i + l) 



£i L £ 2 
L' £i 1 



W 2 )(-l) L V / ^2(4 + l)(2^2 + l) 



^2 L 4 

L' h 1 



<*s,s' , (20) 



{hliSLJW^aiX^ x C (2) ] (1) ||^25'L'J') 

i=l 

= (-1)^-1=^(2 J + i)(2J' + l)(2L + 1)(2L' + 1)(2S + 1)(2S" + 1) | j, f \ j j L L' 2 j 



(«„>( -i.. s " +i '(Mlc- 2 '|K^ * ^ 



W|c (2) ||^) 



K>(-i) s+L (^||c( 2 )||£ 2 ) 





£ 4 1" 




£2 2 / 



/l(l + l)(2l + l) 
(2£ + 3)(2£- 1) 



, (21) 



(22) 



{l 1 t 2 SLJ\\Y j a Si 5 li fiS {1) \\ht2S l L'J') = (-l) s+i+J ' +1 ^ V(2J + 1)(2J' + 1)(2S + 1)(2S' + 1) 



S J L 



S 



ji " s > 1 \\ s> 1 1 M(-l) 5 '(« Sl )^,o + (-l) 6 '(a, 2 )^ 2 ,o]<5L,L' • (23) 
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9.2 Electric quadrupole interaction 



(t&SLJIFlHQie&S'L'J'IF) 



(27 + 3)(2J +!)(/+!) ( J I F 



4/(2/ - 1) 



/ J' 



(^ 2 SL J| |e 2 Q rr 3 c( 2 ) | \l x l 2 S'L'J') , (24) 



(44SL J| |e 2 Q ^ r t - 3 C^ | ^S'L'J') 



i=i 



(- 



. 1)S +L' + J+£ 1+ £ 2 ^ (2J + 1)(2J , + 1)(2L + 1){2U + 1} | £ J 5 j 



(^X-i^Hc' 2 )^) 



L £o 



V 



(b i2 )(-i) L (£ 2 \\c^\\e 2 ) 



£ 2 L t x 
V h 2 



6s,s'- (25) 
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